Tracking Dark Energy with the ISW effect: short and long-term predictions 
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We present an analysis of the constraining power of future measurements of the Integrated Sachs- 
Wolfe (ISW) effect on models of the equation of state of dark energy as a function of redshift, w(z). 
To achieve this, we employ a new parameterization of w, which utilizes the mean value of w(z) 
({w}) as an explicit parameter. This helps to separate the information contained in the estimation 
of the distance to the last scattering surface (from the CMB) from the information contained in 
the ISW effect. We then use Fisher analysis to forecast the expected uncertainties in the measured 
parameters from future ISW observations for two models of dark energy with very different time 
evolution properties. For example, we demonstrate that the cross-correlation of Planck CMB data 
and bSST galaxy catalogs will provide competitive constraints on w(z), compared to a SNAP-like 
SNe project, for models of dark energy with a rapidly changing equation of state (e.g. "Kink" 
models). Our work confirms that, while SNe measurements are more suitable for constraining 
variations in w(z) at low redshift, the ISW effect can provide important independent constraints on 
w(z) at high z. 

PACS numbers: 



I. INTRODUCTION 

There is now a substantial amount of observational 
evidence that the universe is dominated by a dark com- 
ponent which is causing its expansion rate to acceler- 
ate. The analysis of the Cosmic Microwave Background 
(CMB) anisotropy power spectra [l| combined with re- 
sults from large scale structure (LSS) surveys @] strongly 
suggest that about 70% of the energy in the universe is 
in an exotic form of matter which we refer to as dark en- 
ergy (DE). Measurements of the luminosity distance to 
Supernova Type la (SNIa) independently confirm these 
conclusions by showing that high redshift supernovae are 
dimmer than in a matter dominated universe Q . 

Further evidence, which has recently become available, 
is the detection of the Integrated Sachs- Wolfe (ISW) ef- 
fect using the CMB/LSS cross-correlation 0, |E IS 0, 1! 
This evidence is complementary in nature to the SNIa 
data. Rather than probing the overall expansion of the 
universe, it detects the slow down in the growth of den- 
sity perturbations that occurs when the universe ceases 
to be matter dominated. CMB photons traveling to us 
from the surface of last scattering blueshift and redshift 
as they fall into and climb out of gravitational poten- 
tials along their paths. During matter domination, the 
large scale potentials remain constant in time, hence the 
blueshift and the redshift exactly cancel each other out. 
However, any deviation from a constant total equation 
of state results in a time variation of the potentials and 
a net change in the photon energy. This is observed as 
an additional CMB temperature anisotropy, which is the 
ISW effect @. 

Detecting the ISW effect from measurements of CMB 
temperature anisotropy alone is not feasible because the 
signal is hard to separate from the primordial anisotropy 



from the last scattering surface. To circumvent this, it 
was proposed to correlate the CMB temperature maps 
with the local distribution of matter [l(| • The cross- 
correlation arises from CMB anisotropies produced af- 
ter the cosmic structures were formed. On large angular 
scales the cross-correlation signal is dominated by the 
ISW effect, while at small angles an additional contribu- 
tion comes from the Sunyaev-Zel'dovich effect |TTL IT^ . 

The cosmological constant A is the simplest model of 
dark energy that provides a satisfactory fit to the exist- 
ing data. Despite the success of the concordance A cold 
dark matter scenario (ACDM), it is difficult to explain 
why the cosmological constant value is so extremely small 
compared to the expectations of particle physics without 
involving anthropic selection 0] • O n the other hand, the 
observations are also consistent with an evolving dark en- 
ergy, such as Quintessence models where the dark energy 
originates from a scalar field 0, 0, . Establishing 
whether the dark energy is constant or evolving is one 
of the main challenges for modern cosmology. In the 
Friedmann- Robertson- Walker (FRW) universe the evolu- 
tion of dark energy is completely determined by its equa- 
tion of state, which is defined as the ratio of pressure to 
energy density: w = p/ p. For scalar field Quintessence, 
the sound speed is unity and w determines the cluster- 
ing properties of DE. Depending on the model w can be 
constant or change with time. Models with w ^ — 1 cor- 
respond to evolving DE, while w = — 1 corresponds to 
A. 

The time dependence of the dark energy equation has 
been constrained by fitting various forms of w(z) to the 
SNIa data, often in combination with CMB and the LSS 
measurements [H|ll |H E3, El E2 . However the data 
are not accurate enough to distinguish between the cos- 
mological constant and many forms of dynamical dark 
energy. Moreover degeneracies between dark energy pa- 
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rameters strongly limit the possibility to test whether w 
is constant or not. 

As we shall describe in this paper, the CMB/LSS corre- 
lation can provide another probe of the evolution of w(z). 
The sensitivity of the cross-correlation to the dark energy 
evolution has previously been investigated in [23l |24| . In 
this paper we study constraints expected from correlat- 
ing surveys such as DES and LSST with the CMB data. 
We introduce a novel way of parameterizing w(z) that 
has the average (w) (eq. ©) as an explicit parameter. 
This helps to separate information contained in the esti- 
mate of the distance to last scattering from that in the 
ISW contribution. We then proceed to show that the 
cross-correlation, in the absence of errors, is more sensi- 
tive to the details of the high redshift evolution of w(z) 
than other available observables. This sensitivity, how- 
ever, is off-set by large statistical errors due to a large 
primordial contribution to the CMB signal. We forecast 
the expected uncertainties in DE parameters (using two 
different models) extracted from CMB/LSS correlation 
that will be possible in the next five years (WMAP/DES) 
and in the next ten years (Planck/LSST). We also calcu- 
late corresponding constraints from the SNe data (SNLS 
in short term and SNAP in long term). We find that 
in long term, cross-correlation can provide competitive 
constraints on the evolution of w(z). 

The paper is organized as follows. In Section [n] we 
discuss the dark energy modeling and the details of the 
calculations. Section UTTl describes the observations con- 
sidered in the paper. In Section Hvl we discuss some ad- 
vantages of including the cross correlation information 
and give the resulting parameter constraints. Finally in 
Section we present our conclusions. 



II. FORMALISM AND METHODS 
A. Dark energy model 

Many models of dark energy have been proposed: cos- 
mological constant (A), Quintessence, k-essence, tangled 
defects, etc. For simplicity, we will focus on quintessence 
based models where the dark energy arises from a scalar 
field currently slow-rolling down a potential, which is ef- 
fectively inflation occurring at late times. Perturbations 
in the dark energy are easily incorporated in such mod- 
els, and the sound speed is typically of order the speed 
of light. 

The shape of the potential, which determines the evo- 
lution of the dark energy, is generally a free function. 
This makes it difficult to predict the time dependence of 
the DE equation of state, so one is potentially left with 
constraining a completely arbitrary function of redshift. 
The problem must be simplified by considering some pa- 
rameterization of the function w(z). Several parameteri- 
zations have been suggested in a vast literature, but any 
resulting constraints are very sensitive to the choice of 
parameterization. For instance it was shown in |25j | that 



fitting a constant w will tend to favor lower negative val- 
ues if the DE equation of state is time dependent. Even 
two-parameter Taylor expansions are too simplistic, re- 
sulting in a strong dependence on the assum ed p riors and 
generally biased against dynamical models [171 l20l |26( . 

One of the most popular two-parameter formulae is the 
linear change in the scale factor a = (1 + z)^ 1 [27j given 
by, 



w(a) = wq + (a — l)S n 



(1) 



where wq = wq (today) and w rn = wq + S w being the 
value at some early time, such as the radiation-matter 
equality. From eq. it appears evident that with 
two parameters one can fix the value of the equation 
of state today and at high redshift, but the time evo- 
lution between the two extremes cannot be changed (i.e. 
dw/da — const.). As consequence of this, a linear ex- 
pansion cannot allow for rapid variation in the equation 
of state and particularly in the case of Quintessence, for 
which wq > — 1, this bias will tend to favor small values 
of S w . Despite these potential problems, eq. Q has be- 
come a common tool of testing the capability of future 
experiments to constrain DE models (see |28|). Due to 
its wide use in the literature, we will consider this pa- 
rameterization here to help the comparison of our results 
with previous works. We shall refer to it as Model I. 

In order to avoid the pitfalls introduced by the Taylor 
expansion we also consider a phenomenological approach 
and model w(z) according to some minimal requirements, 
which are: 1) the form of w{z) depends on a minimal 
number of parameters which are of simple physical inter- 
pretation, 2) the parameterization is not biased against 
any particular time behavior and can account for both 
rapid (large) or slow (small) variation of the equation of 
state, 3) the parameterization reproduces the evolution 
of several proposed models of DE such as Quintessence. 
In this regard, a general formula constructed to follow 
these guidelines was proposed in I29I l30| and is usually 
referred to as the "Kink" model [lg . 

Within the class of the Kink models we consider a 
unique new parameterization. We start with the func- 
tional form 



w{z) 



W + + WJ_ w + 



tanh 



Z T 



(2) 



which describes a transition from w+ for z — ► +oo, to 
w_ for z — > — oo, with parameters A z and zt describing 
the width and the central redshift of the transition. The 
novelty is in using the average equation of state ((w)) as 
an explicit parameter, where (w) is defined as 



f 1 da n D (a) 



(3) 



Here a is the scale factor (a = 1 today), a^s is the scale 
factor at the surface of last scattering and fin (a) is the 
ratio of the DE energy density to the critical density. 
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FIG. 1: The equation of state as a function of redshift for 
the two fiducial models considered in the paper. Model I 
(blue solid line) and Model II (black dash line) have the same 
weighted average (w) = —0.9 (the dotted line) defined in 
eq. Q. 



It is well known [3ll 132 ] that on small scales the CMB 
temperature anisotropy spectra are insensitive to the de- 
tails of w(z) except for its average value as defined in 
eq. J3J). Specifying (w) roughly fixes the distance to the 
last scattering surface |33| , which determines positions of 
the acoustic peaks in the CMB temperature anisotropy 
spectrum. The CMB temperature data, combined with 
a prior knowledge of Hq and under the assumption that 
the universe is flat, is capable of constraining the value of 
(w) with a very high accuracy. By explicitly specifying 
(w) we preserve the information that would be obtained 
from fitting a constant w to the CMB, while still allowing 
for a varying w(z). Hence, a direct control on the value 
(w) allows us to reduce the four dimensional parameter 
space to a much smaller observationally allowed region, 
which corresponds to a narrow range of possible values 
of (w) . 

Note that for (w) close to —1, as favored by the data, 
most choices of A. w would lead to w{z) changing between 
— 1 and — 1 + 1 A. w | , where A^, = u>_ — u>+ . Also note that 
in a large class of quintessence models either w + or w- 
is —1. This is the case when the scalar field is initially in 
a slow-roll reg ime {w+ = —1) as in, e.g., the Doomsday 
model |33. l35l|. Similarly for tracker models, where the 
quintessence field starts tracking the background com- 
ponent and evolves at late time settling to a minimum 
of its potential (i. e. w- = — 1) 36]. Hence, fixing the 
value of either w_ or u> + to —1 is physically and obser- 
vationally motivated and allows us to reduce the number 
of parameters in eq. J2J from four to three. 

Thus, instead of using the four parameters of eq. 



w+ , W-, zt and A z , we use three parameters: (w), A w 
and A z . Here, zt is chosen to reproduce the desired 
value of (w). (Alternatively, we could explicitly specify 
zt and use eq. © to find the corresponding value of A z .) 
Since we are primarily interested in determining whether 
w(z) is varying or constant, A w is the main parameter 
of interest. The role of the third parameter, zt or A z , is 
only to allow for sufficient freedom in the ansatz. In this 
paper, we will stay with the choice of (w), A w and A z 
and refer to this parameterization as our Model II. 

Representative plots of w(z) versus z for Model I and 
Model II are shown in Fig. ^ 



B. Temperature and density correlations 

The CMB temperature anisotropy due to the ISW ef- 
fect can be written as 



P° ( ■ 

\ drj ( $ [(ry - r?)n, 77] - * [(r; - 77)11, 77] 



ST(A) 
T 

.(4) 

where n is a direction on the sky, 77 is the conformal time, 
77 m is some initial time far into the matter era, 770 is the 
time today, $ and '5 are the gravitational potentials in 
the Newtonian gauge, and the dot denotes differentiation 
with respect to 77. 

The ISW temperature anisotropy can be correlated 
with the distribution of galaxies on the sky using, 



Sg(p) 



2V(n) - N 
N 



(5) 



where S g (fi) is the overdensity of galaxies in the direc- 
tion n, N(n) is the number of galaxies in the pixel corre- 
sponding to the direction n, and N is the mean number 
of galaxies per pixel. We assume here that we sample 
galaxies in a fixed redshift region, characterized by a nor- 
malized galaxy selection function, W g (z). For simplicity, 
we initially consider a single selection function, but below 
we consider correlations arising between different redshift 
bins as well. 

The galaxy number overdensity S g (h) is assumed to 
be tracing the cold dark matter (CDM) density contrast 
<5 c (n) via 



b g S c (n), 



(6) 



where b g is the linear galaxy bias, which is possibly red- 
shift dependent. The CDM density contrast can be writ- 
ten as an integral over conformal time, 



£ c (n) 



n ° dz 



77)11,77), (7) 



where (5 c (x,77) is the three-dimensional density contrast 
and is directly related to the gravitational potentials. 

We are interested in calculating the cross-correlation 
function 



/ ST 1SW 

X{9) = XQ&i n 2 |) = ^__(n l)(5s (n 2 ) 



(8) 
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where the angular brackets denote ensemble averaging 
and 9 is the angle between directions ni and n 2 . Cal- 
culations are often simplified when X(9) is decomposed 
into a Legendre series, 

X{6) - Y, ^ ± X l P l {cos 9) . (9) 

1=2 7F 

The coefficients Xg can be evaluated using j^: 

Xl = 4n ~kjt An{k) A ^ W(fc) Me{k) ' (10) 

where k is the wave- number, A^(fc) is the primordial 
curvature power spectrum, and functions A* sw (fc) and 
Mi(k) are defined as 

Af w (fc) = / d V j e (k[v-Vo}) 
x {c^ys(j){k,r}) - ip(k,r))) 

l-rio 

Mt(k) = bgCgy / dr) ji(k[r] - rjo]) 

x zW g {z{ri))8{k,ri) . (ff) 

In the above, the integration starts at an early time 
when a given mode k is well-outside the horizon, ji(-) are 
the spherical Bessel functions, </>(fc, ry), ^(k, rj) and S(k, rf) 
are the evolution functions defined in |23j along with the 
numerical coefficients and c$^,: 

6 3 _ $ / 2 \ 

c<5 * = vj/ = ~2 ' C ** = * = ~ \ 1 + 5 j ' ^ ' 

where i?„ = p„/{p 1 + p v ) and p 1 and are the photon 
and relativistic neutrino densities. 

Above we have considered the correlations between the 
CMB and a sample of galaxies defined by a single galaxy 
selection function. Future surveys will enable us to sep- 
arate galaxies into many redshift slices with photometric 
redshift slices (bins), and we can consider separate corre- 
lations with each of the slices, X\. We calculate these as 
above using the different selection functions, Wg{z), each 
of which has a corresponding weighting function Mj(k). 
For each bin, we consider a different possible bias fac- 
tor b g , but as we discuss in Section HV Bl these can be 
determined sufficiently well by the observations of the 
galaxy-galaxy correlation functions. 

The galaxy-galaxy correlations are evaluated in a sim- 
ilar way. We have many correlation functions w^'(6), 
corresponding to correlations between all possible red- 
shift bins labeled by indices Similarly to eq. QlOjl. 
the Legendre coefficients of these correlations are given 
by 

The full CMB auto-correlation Cj T contains the pri- 
mary anisotropy from the last scattering surface as well 



as contributions from the ISW, the SZ and other effects: 

Cj T = Cf s + Cf w + Cf z + ... . (14) 

The uncorrelated anisotropics from the last scattering 
surface dominate the noise in the cross correlation de- 
tection. The ISW contribution dominates the galaxy- 
CMB cross-correlation on large scales (where the signal 
is strongest), while the SZ effect is anti-correlated at 
WMAP frequencies and affects only the small angular 
scales [6|- In principle the SZ signal can be eliminated 
by smoothing the CMB maps on scales smaller than the 
typical angular size of a cluster ( 1°), therefore we will 
not include the SZ in our analysis. 

We evaluate the CMB/LSS correlation coefficients Xj, 
as well as all other relevant spectra, using an appropri- 
ately modified version of CMBFAST 37J. For all models, 
we assume a flat universe with adiabatic initial conditions 
for all particle species, including the Quintessence. 



C. Fisher matrices 

We use the usual Fisher method [3^| for parameter esti- 
mation forecasts to study the potential of the CMB/LSS 
cross-correlation for dark energy constraints. We as- 
sume that the universe has evolved from adiabatic ini- 
tial conditions with a nearly scale-invariant spectrum of 
density fluctuation, negligible tensor component and no 
massive neutrinos. It is also assumed to be flat with 
a Quintessence field parameterized by w{z). In addi- 
tion to the dark energy parameters [wo^to] of Model I 
or [(w) ,A W ,A Z ] of Model II, our parameter space p a in- 
cludes the matter fraction of total energy density Hm, 
the Hubble parameter h, the baryon density u>b — Qbh 2 , 
the reionization optical depth r, the scalar spectral in- 
dex n s and the amplitude of the initial power spectrum 
A s . In total, we have 8 parameters for Model I and 9 pa- 
rameters for Model II. In addition, for the SNe analysis, 
we account for M - the intrinsic supernova magnitude. 
Cross-correlation parameters also include the bias factors 
for each photometric bin, as described in Section TlV Bl 

The cross-correlation Fisher matrix can be written as 

^ = ^ ( £EScov-.(XlX»g. ,15) 

Here Cov~ 1 (X]AT^), as well as the partial derivatives, are 
evaluated at the fiducial values (see Sec. Ill L) I) . The sky 
fraction, / s k y , is the smallest of the CMB and the galaxy 
survey sky coverage fractions. We used £x = 1000 as 
the summation limit, although most of the contribution 
comes from £ < 500. The smallest value of I can be 
approximately determined from £ m ; n « 7r/(2/ s k y ). The 
covariance matrix Cov^ is given by 

Cov(XjX|) = (cf T 4 i,j) + X}Xj) j {21 + 1) , (16) 
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where Cj T and are denned in eqns. (|3"0|) and if2*5|) 

and include the contributions from the noise. Note that 
the contribution to the Cj T from the surface of last scat- 
tering is much larger than that from the ISW effect, while 
Xf is not sensitive to the last scattering physics at all. 
This results in a very large variance in Xj that signifi- 
cantly impairs its potential for parameter estimation. 

For CMB measurements, the Fisher information ma- 
trix is given by 



^CMB 

*a,0 ~ /sky 2^ 2^ :% 
1=2 A,B UP ' 



Cov-^Cf)^, (17) 



•a 



with the covariance matrices for A,B = TT, EE, TE 
given explicitly in 39]. In the above, ^cmb was taken 
to be 1200 for WMAP and 2000 for Planck, safely above 
the maximum I allowed by the angular resolution of each 
experiment. The partial derivatives and the covariance 
matrix are evaluated at a fiducial choice of parameters 
that is specified in the next subsection. 

The quantity directly measured from SNe observations 
is their redshift-dependent magnitude 




m(z) = M + b\ogd L + 25 



(18) 



FIG. 2: CMB temperature angular power spectrum (TT) and 
the temperature-polarization cross-correlation (TE) for the 
three models in Fig. Q (using the same conventions for the 
three line types) with the WMAP's first year data. 



where M is the intrinsic supernova magnitude and dh the 
luminosity distance (in Mpc). The luminosity distance is 
defined as 



d L {z) = {\ + z) 



dz' 



(19) 



where H(z) is the Hubble parameter with a current value 
of Hq. The information matrix for SNe observations is 



JV 

P SN _ \ " 



1 dm(zi) dm(zi) 
a(zi) 2 dp a dp/3 



(20) 



where the summation is over the redshift bins and a m [zi) 
is the value given by eq. (|31[) at the midpoint of the i-th 
bin. 

Given the Fisher matrix, the lower bound on the un- 
certainty in determination of a parameter pi is given by 



(21) 




k (h/Mpc) 



where Cy is the error matrix. 

One can forecast constraints from a combination of 
measurements by adding the individual Fisher matri- 
ces. For example, cross-correlation measurements can 
be combined with measurements of the CMB spectra by 
Planck as 



rpPlanck+X 

a, (3 



l^Planck 



■ F a,f3 



(22) 



D. The fiducial models 



Forecasts obtained using Fisher analysis are sensitive 
to the choice of the fiducial model. We will consider the 



FIG. 3: The linear matter power spectrum for the three mod- 
els in Fig. The SDSS data points |4l| are plotted to show 
current experimental error bars. 



same fiducial models for the two dark energy parameter- 
izations discussed above and plotted in Fig. ^ The fidu- 
cial parameters for Model I were taken to be wq = —0.98 
and 8 W = -0.29. For Model II we chose (w) = -0.9, 
A w = —0.6 and A 2 = 0.3. We take all other cosmo- 
logical parameters to be the same for both models. We 
have set the Hubble parameter h = 0.69, baryon density 
fl b h 2 = 0.024, total matter density fl M h 2 = 0.14, 



spec- 
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FIG. 4: The percent difference in Ho<1l{z) between the model 
with constant w = —0.9 (dotted line) and Models I (solid 
blue) and Model II (black dash). Current SNa data can de- 
termine HocLl(z) with an accuracy of 5 — 10%. Future data 
sets, such as supernovae from SNAP, can achieve an accuracy 
of 1 - 2%. 



tral index n s — 0.99, optical depth t = 0.166 and the 
amplitude of scalar fluctuations A s = 0.86 (as defined in 

m) 

All of current data is consistent with w = const = —0.9 
at la level Q. Both, Model I and Model II, have the 
same weighted average (w) = —0.9. This assures that 
these two models have nearly identical CMB spectra as 
shown in Fig. |2J 

In Fig. |3 we compare the matter power spectra pre- 
dicted by the two models. Plotted are the linear CDM 
spectra at z = 0. The bias, non-linear effects and the 
redshift space distortion would modify the three spectra 
in a similar way. While one cannot compare these linear 
spectra to the data without accounting for the above- 
mentioned corrections, we include the SDSS data points 
|4lj to illustrate the uncertainty in the current determi- 
nation of P(k). 

Models I and II are also consistent with current es- 
timates of the change in luminosity distance d^iz) ob- 
tained from the SNIa measurements. In Fig. 21 we plot 
the percent difference in HqcIl between the model with 
constant w — —0.9 and Models I and Model II. The dif- 
ferences between the models are of order 1 — 2%, well 
below the current level uncertainty and comparable to 
the projected accuracy of SNAP. 

Hence, present data cannot distinguish between w = 
const = —0.9 and our Models I and II and both are 
perfectly consistent with all available observations. 
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FIG. 5: The galaxy number distribution vs. redshift expected 
from the DES survey. 



III. EXPERIMENTS 

We are interested in forecasting the errors in dark en- 
ergy parameters that can be extracted from CMB/LSS 
correlation studies and compare them to those expected 
from the luminosity distance measurements. We make 
short- (less than 5 years) and long- (ten years) term pre- 
dictions based on the ongoing and planned CMB, galaxy 
and SNe observations. 

For our short term predictions we assume the 4-year 
CMB temperature and polarization data from WMAP 
|42|, the expected galaxy counts from the Dark Energy 
Survey (DES) [13], and the SNe from the Nearby Super- 
nova Factory (NSNF) |44j and the Supernovae Legacy 
Survey (SNLS) [45|. In addition, we impose a Gaussian 
prior on the value of h corresponding to a la uncertainty 
of ±0.08 from the Hubble Space Telescope's Key Project 
(HSTKP) HI. 

For the long term predictions we assume the 14 months 
CMB temperature and polarization data by Planck satel- 
lite 01 1 the galaxy catalogues by the Large Synoptic 
Survey Telescope (LSST) [13], and the SNe from the Su- 
pernovae Acceleration Probe (SNAP) |49j complemented 
by those from the NSNF. 



A. Galaxy data: DES and LSST 

To describe the distribution of galaxies from a given 
survey, we take the total galaxy number density to have 
the following redshift dependence [Hcj : 

n t g ot (z) C xz 2 e-^ z ^\ (23) 
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tometric redshift distributions are given by, 



1 



<(z) = -n\°\z) 



erfc 



Zi-l 



( V- 



erfc 



/ Zj - z 
\V2a(z) 



(25) 

where erfc is the complementary error function. 

For a given photometric redshift bin, the normalized 
selection function is given by 



wi{z) 



n l(z) 



(26) 



where N l is the total number of galaxies in the i-th bin. 
The Poisson noise is uncorrelated between bins and con- 
tributes only to the galaxy angular auto-correlation as 



N (id) = £y 



(27) 



FIG. 6: The galaxy number distribution vs. redshift for the 
"conservative" LSST survey. 



where Sij is the Kronecker delta and n l A is the galaxy 
number per solid angle in the i-th bin. The observed 
correlation u>) is the sum of the signal and the Poisson 




Aid) 



(id) 



(28) 



FIG. 7: The galaxy number distribution vs. redshift for the 
"goal" LSST survey. 



where z n is a parameter that is adjusted to reproduce the 
expected median redshift of the survey. The galaxies are 
further divided into photometric redshift bins, i. e. 



'(*)=XX(*)- 



(24) 



Following reference [50|. we assume that the photomet- 
ric redshift errors are Gaussian distributed, and that 
their rms fluctuations increase with redshift as <r{z) = 
c max (l + z)/(l + z m ax)- The bins sizes are chosen to 
increase proportionally to the errors. The resulting pho- 



The Dark Energy Survey (DES) is designed to probe 
the redshift range 0.1 < z < 1.3 with an approximate 
l-er error of 0.1 in photometric redshift. This roughly 
corresponds to four bins (see Fig. |3J and a median red- 
shift of z = 0.6. The total expected number of galaxies 
is approximately 250 million in a 5000 sq. deg. area on 
the sky, or f sky = 0.13. 

The proposed LSST survey is expected to cover up to 
20, 000 sq. deg. of the sky and catalogue over a billion 
galaxies out to z ~ 3. We consider two cases: the conser- 
vative scenario case with f s k y = 0.3 and 7 photometric 
redshift bins out to z ~ 1.5, and the desired case with 
fsky = 0.5 and 10 photometric redshift bins out to z ~ 3. 
For both cases we assume 70 gal/arcmin 2 . The assumed 
galaxy number distributions for the two cases are shown 
in Figs. E]and[7| In what follows, we will refer to the two 
LSST cases as 11 conservative" and "goaF . 

The photometric redshift errors used in our analysis 
should be seen as the optimistic values for the respective 
experiments. For a discussion of potential sources of sys- 
tematic errors in photo-z estimates the reader is referred 

to mm. 



B. CMB data: WMAP and Planck 

The latest expected values for the sensitivity and 
the resolution parameters of the 4-year WMAP and 14 
months Planck missions are available in E3- To 
make comparison with previous results easier, we opted 
to use similar, although somewhat different, parameters 
employed in |53ll54|| for the three highest frequency chan- 
nels of WMAP and the lowest three Planck HFI channels. 
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We have checked that using the exact parameters given 
in 03 would lead to an insignificant modification of 
our results. The relevant parameters are listed in Ta- 
ble [I] The noise contribution to the CMB temperature 

TABLE I: The relevant parameters of WMAP g3 and Planck 
jr^l . We use the three highest frequency WMAP channels and 
the lowest Planck HFI frequency channels. 



WMAP Planck 



v (GHz) 


41 


61 


94 


100 


143 


217 


Afwhm (arc min) 


28 


21 


13 


10.7 


8.0 


5.5 


<JT (/iK) 


22 


29 


49 


5.4 


6.0 


13.1 


<TE (a*K) 


30 


45 


75 


n/a 


11.4 


26.7 


/sky 




0.8 






0.8 





The SNLS is expected to measure approximately 700 
supernovae out to z max = 1 with an uncertainty 5m — 
0.02. The SNLS supernovae distribution that we use 
for our short term forecasts is shown in Tabic [H] The 
ground-based NSNF experiment |44| will add to the 
count 300 SNe at z < 0.1. 

SNAP will provide over 2000 supernovae with 8m = 
0.02 and z max = 1.7. In this work we have considered a 
fiducial SNAP survey analogous to the one modeled in 
The parameters are given in Table ITT1 

In addition to supernovae, SNAP can get over 300 mil- 
lion galaxies (at n ga i = 100 arcmin -2 ) at moderate and 
high redshifts, covering potentially up to 10, 000 sq.deg 
of the sky. This would make it a candidate survey for 
ISW studies 1 . 



auto-correlation (TT) and the E-mode polarization (EE) 
spectra from one frequency channel is given by 



N, 



AA 



VA,c PFWHM, 
TcMB 



„/81n2 



where c labels the channel and A — T,E. The combined 
noise from all channels is 



N AA = 



(29) 



The observed spectra, i.e. the signal plus the noise, are 
then 



Cf A = Cf A + Nf A . 



SN data: SNLS and SNAP 



(30) 



The DE constraints from supernovae surveys depend 
on the depth of the survey, the total number of the super- 
novae and their redshift distribution, i. e. on the number 
of SNe per redshift bin. For example, when constrain- 
ing an evolving equation of state, it is known |55| that 
the choice of the distribution affects the bounds on the 
relevant parameters (e.g. [wq,S w ]). 

SNe observations determine the magnitude, m(z), de- 
fined in eq. (|18l) . The uncertainty in m(z), at any z-bin 
containing A^m supernovae, is given by 



'obs 



N h 



dm 2 , 



(31) 



and we assume <7 bs = 0.15. The systematic error, dm, 
is assumed to increase linearly with redshift: 



dm = 5r 



(32) 



5m being the expected uncertainty and z max the maxi- 
mum redshift. 



IV. RESULTS AND DISCUSSION 
A. Cross correlations and w(z) 

As discussed earlier, CMB is sensitive primarily to a 
particular weighted avera ge o f the dark energy equation 
of state, given by (w) [3l|, |32, |33 ■ Measurements of the 
acceleration by SNe are predominantly sensitive to the 
very low redshift behavior of w(z). The cross correlation 
measurements, on the other hand, can give an indepen- 
dent window on the behavior of the equation of state at 
intermediate redshifts. In this subsection we justify our 
expectation for the CMB/LSS cross-correlation to be an 
effective probe of w(z). Additional analysis of the de- 
pendence of the cross correlation on the time evolution 
of the DE equation of state can be found in [2^, 01 • 

First note that in the small angle limit, the cross cor- 
relation can be written as 



X 



(0 -> 0) rs const [ dz WJz)D(z) — [(1 + z)D(z)] 
J dz 



(33) 

where D(z) = S(k, z)/6(k, 0) is the linear growth factor 
and we used the Poisson equation to express $ in terms 
of 5: &(k,z) oc (1 + z)S(k,z). The small angle limit 
X(9 — > 0) is representative of the total cross-correlation 
signal inside the redshift window W g (z) 24]. 

Eq. (|33|l shows that the cross-correlation is essentially 
the product of the growth function and its derivative av- 
eraged over a given range of redshifts. Having several 
weakly overlapping selection functions Wt(z), as in the 
case of LSST (e.g. see Fig. 0, would allow one to effec- 
tively map the evolution of this product. Since the rate of 
the growth is a particularly sensitive probe of the details 
of the dark energy evolution (see [53 for a discussion), 
cross-correlation studies can provide useful information 
despite large statistical errors. To illustrate the sensitiv- 



1 We thank Eric Linder for pointing this out 
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TABLE II: The redshift distribution of type la supernovae N(z) for SNLS ^ and SNAP together with 300 SNe from the 
NSNF. The redshifts given are the upper limits of each bin. Magnitude errors cr m (z) are evaluated at bin midpoints. 



z — > 


0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1.0 1.1 1.2 1.3 1.4 1.5 1.6 1.7 


SNLS 


N(z) 

a m (z) (xlO" 3 ) 


300 56 70 84 133 105 84 84 42 7 
9 20 19 18 16 19 22 23 29 60 


SNAP 


N(z) 

a m {z) (xKT 3 ) 


300 35 64 95 124 150 171 183 179 170 155 142 130 119 107 94 80 
9 25 19 16 14 14 14 14 15 16 17 18 20 21 22 23 26 




where w(z) is not affected by a change in A w . That is, 
w = —1 at low redshifts for any small variation around 
Model II. 

A change in (while holding (w) fixed) affects the 
observables via a combination of two effects. One is the 
change in the high redshift value of w, which alters the 
DE fraction at early times. The second is the change 
in the transition redshift, zt- The change to the growth 
factor comes mainly from the first of the two effects, while 
P{z) is affected by both. 

The sensitivity of P{z) to w{z) is off-set by the large 
variance in cross-correlation measurements. Below we 
present results of the Fisher analysis that reflect this lim- 
itation. 



j i i i i i i i i i i 

12 3 



FIG. 8: The percent change in P(z) (eq. 1341 1 . solid line), 
dh{z) (dot), [D(z)] 2 (dash) and SIde(z) (dash-dot) caused by 
a 10% change in A w , while holding (w) fixed. 

ity of the cross-correlation to changes in the evolution of 
w(z), we look at the change in the quantity 

P(z) = D(z)^[(l + z)D{z)} (34) 

caused by a 10% decrease in the parameter A w of our 
Model II. We hold the other two DE parameters, (w) and 
A z , as well as the cosmological parameters, fixed at their 
fiducial values. Then we compare the resulting relative 
difference in P(z) to the differences in \D(z)] 2 and in 
cIl(z). As one can see in Fig. [SJ P(z) is by far the more 
sensitive probe of the three. Also note, that the evolution 
of P(z) is almost identically tracing the evolution of the 
DE fraction Qde(z)- This is expected, since the change 
in the gravitational potential $(z) is directly caused by 
a non-zero Qde(z). In particular, it is shown in |Hsj that 
the evolution of (1 + z)D(z) couples to w(z)Q,de{z). 

It is possible to explain the differences in Fig. [S] qual- 
itatively. The reason cIl(z) is not sensitive to a change 
in A w is mainly due to the change occurring at a high 
redshift: zj* ~ 1.2. The main contribution to <Il{z) in 
the integral given by eq. I|19[) comes from low redshifts, 



B. Bias 

Most of the contribution to the CMB/LSS correlation 
signal comes from large scales 2 , where perturbations are 
well described by the linear theory. On such scales, galax- 
ies are expected to closely trace the distribution of dark 
matter, up to a scale- independent bias factor b g . The 
bias is also known to vary with redshift, hence, each of 
the photometric bins will, in principle, have a different 
bias factor. The bias factors corresponding to each bin 
can be determined from the amplitude of the primordial 
spectrum extracted from the CMB and the galaxy-galaxy 
autocorrelation spectra. It is reasonable to assume that 
the bias in each bin can be determined with 10% accu- 
racy. 

We account for the bias uncertainty by assigning an in- 
dependent constant bias factor to each bin and treating 
them as parameters in our Fisher analysis. We assume 
that the value of the bias can be determined from else- 
where to a 10% accuracy in all bins. This assumption, 



2 The actual physical scale depends on the redshift of the selec- 
tion function and on the model. For reference, for galaxies at 
z ~ 0.2 the dominant cross-correlation s cale for the model with 
a constant w = —0.9 is around 50Mpc/h l23l . On smaller scales 
the ISW effect is negligible (the potential wells have to be long 
enough for photons to notice the effect of stretching by the ac- 
celerated expansion). On very large scales the cross-correlation 
signal naturally decreases with the reduction in clustering. 
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FIG. 9: The dependence of the uncertainty in the Model II 
parameter A m on the assumed bias prior. The long dash line 
corresponds to the error from the CMB/LSS cross-correlation 
using WMAP and DES, the short dash-line from Planck and 
the "conservative" LSST, and the solid line from Planck and 
the "goal" LSST. The blue dot lines are the expected errors in 
A™ from the WMAP data alone (the upper dotted line) and 
from Planck alone (the lower dotted line) shown for reference. 

which wc implement by imposing a prior, does not have 
a big influence on our final results, as long as the bias in 
each bin is known to better than 20%. This is illustrated 
in Fig. [5J where we plot the predicted Fisher error in A w 
as a function of the assumed uncertainty in bias. The 
fiducial values of the bias parameters were chosen to be 
the same for all bins. Choosing them according to a fixed 
function of redshift does not noticeably affect the results 
of the Fisher analysis. 

C. Short term prospects 

In the next few years we will have results from DES, 
SNLS and the complete data from WMAP, which will 
include E-modc polarization. In this subsection we show 
and compare the expected constraints on w(z) from these 
measurements . 

The linear parameterization of w(z) (eq. (JQ) with 
fiducial parameters corresponding to ACDM (wq = — 1, 
S-w = 0) is the most commonly considered case in the 
existing literature. For that reason we include it in our 
analysis and show the short term prediction for the DE 
parameters in Fig. 1101 Shown are only the contours for 
Qm> w q an d S w , with other parameters being marginal- 
ized over. In this case, the cross-correlation informa- 
tion provides only a modest improvement over the CMB 
spectra alone. Supernovae data, on the other hand, can 
tighten the allowed range of the DE model parameters 




-2 -1.5 -1 -0.5 




-2 -1.5 -1 -0.5 -2 -10 12 



FIG. 10: The expected la projected contours for the ACDM 
fiducial model from WMAP (thin blue), and the informa- 
tion from WMAP combined with the WMAP/DES correla- 
tion (thick solid red) and SNLS supernovae (black dot), using 
the linear parameterization of eq. (0. For each line-type, the 
smaller ellipses correspond to including the CMB polarization 
information, while the larger ellipses are obtained using the 
CMB temperature spectra only. A prior from HSTKP on the 
value of h is used. 
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FIG. 11: Short term forecasts of uncertainties on the param- 
eters of Model I. Shown are the la projected contours from 
WMAP (thin blue), and the WMAP information combined 
with the WMAP/DES correlation (thick solid red) and SNLS 
supernovae (black dot). As in Fig. 1101 the smaller ellipses 
correspond to including the CMB polarization information, 
while the larger ellipses are obtained using the CMB temper- 
ature spectra only. A prior from HSTKP on the value of h is 
used. 
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FIG. 12: The expected la projected contours from WMAP 
(thin blue), and the WMAP information combined with the 
WMAP/DES correlation (thick solid red) and SNLS super- 
novae (black dot) for Model II. The smaller ellipses corre- 
spond to including the CMB polarization information, while 
the larger ellipses are obtained from the CMB temperature 
spectra alone. A prior from HSTKP on the value of h is used. 





0.295 
0.29 
0.285 



FIG. 13: The expected la projected contours from Planck 
(thin blue), and Planck information combined with the 
Planck/LSST(goal) correlation (thick solid red) and SNAP 
supernovae (black dot) for the linear parameterization of 
eq. Q with ACDM fiducial parameters. (Compare this to 
Fig.Hnj. 



somewhat better and significantly constrain wq. 

The short term predictions for Model I are shown in 
Fig. ^2 again, for Qm, and 5 W only. The list of 
predicted uncertainties in all parameters of Model I is 
given in Table II I II The expected uncertainties are dif- 
ferent, and generally larger than in the ACDM fiducial 
case. This illustrates the importance of the underlying 
DE model for constraining w{z). The cross-correlation 
will not, in short term, provide competitive constraints 
on the evolution of w{z) according to ACDM or Model I. 

The forecasts for Model II parameters (w) , A w and FIm 
are shown in Fig.^J while Table ITTT1 contains the full list 
of the expected parameter uncertainties. The choice of 
the model clearly makes a big difference. Most notably, 
neither the short term SNe nor the cross-correlation mea- 
surements can improve the CMB bounds on A w . The 
SNe, on the other hand, are still useful in tightening 
the constraint on Am and (w). Overall, it is clear from 
the plots and the table that measurements of the cross- 
correlation in the next five years will not provide any 
new competitive constraint on the evolution of w. In the 
case of Model II this is primarily due to limited depth 
of the DES survey, which does not sample the redshifts 
above z ~ 1, where the transition in w(z) occurs. In the 
case of Model I, constraints from the cross-correlation 
are even weaker, even though it has one less parameter 
as compared to Model II. It is worth noticing that for 
this model, the degeneracy between wq and 5 W that hin- 
ders their determination from CMB alone (with a prior 
on h) is slightly improved by the addition of the cross- 
correlation information. This is because the angular di- 
ameter distance to last scattering surface depends on the 
dark energy through the average equation of state value 



(w) . Therefore in Model I different values of wq and 5 W 
can give the same (w). On the other hand redshift depen- 
dent measurements of the ISW-correlation are sensitive 
to both the average equation of state and its time evolu- 
tion improving the constraints on these two parameters. 

In summary, the short term cross-correlation measure- 
ments will not be competitive, while supernovae mea- 
surements will only constrain wq of Model I and (w) of 
Model II. As we will see in the next subsection, the con- 
straints on all parameters tighten in the case of deeper 
surveys such as LSST. 



D. Long term prospects 

In the longer term, i.e. in the next ten years, we should 
expect to have the results from Planck CMB measure- 
ments, galaxy catalogues from LSST and the SNe mea- 
surements from SNAP. These will significantly improve 
the constraints on the DE parameters. In certain mod- 
els, as we show here on the example of our Model II, 
CMB/LSS cross-correlation can play a competitive role. 

As in the previous subsection, we begin by considering 
the linear DE model with two different sets of the fiducial 
values: one where they match a cosmological constant 
(wq = — 1, 6 W = 0) and another where w evolves but still 
consistent with present observations (Model I). In Fig. 1131 
and Fig. ^] we show the la projected contours for the 
linear models. On these plots, the cross-correlation and 
the SNe prediction includes the information from Planck. 
Combining with Planck is equivalent to applying strong 
priors on all cosmological parameters. The constraints 
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TABLE III: Short term predictions assuming WMAP for CMB, SNLS+NSNF for SNe and WMAP/DES for the CMB/LSS 
cross- correlation (X). The SNe and X results are shown with priors from CMB temperature only (T) and CMB temperature 
and polarization data combined (TP). A prior from HSTKP on the value of h is used. 







Model I, la errors 


Model II, la errors 






CMB 


SNe 


X 


CMB 


SNe 


X 


p 


fiducial 

ValllCR 


T 


TP 


T 


TP 


T 


TP 


T 


TP 


T 


TP 


T 


TP 


w 


-0.98 


1.1 


0.73 


0.16 


0.15 


0.87 


0.62 


~ 


~ 


~ 


- 


~ 


- 


5 W 


-0.29 


3.8 


2.1 


0.82 


0.66 


3.0 


1.8 














(w) 


-0.9 














0.26 


0.19 


0.22 


0.13 


0.23 


0.17 


A w 


-0.6 














3.2 


1.0 


2.9 


1.0 


2.0 


0.97 


A z 


0.3 














8.9 


6.2 


4.0 


2.6 


6.6 


4.7 




0.3 


0.07 


0.06 


0.04 


0.02 


0.07 


0.06 


0.07 


0.06 


0.05 


0.03 


0.07 


0.05 


h 


0.69 


0.07 


0.07 


0.03 


0.02 


0.07 


0.07 


0.08 


0.07 


0.06 


0.04 


0.07 


0.06 


10 3 uj b 


24 


2 


0.7 


2 


0.6 


2 


0.7 


3.4 


0.8 


3.2 


0.7 


2.7 


0.7 


lis 


0.99 


0.08 


0.02 


0.07 


0.02 


0.07 


0.02 


0.13 


0.02 


0.12 


0.02 


0.1 


0.02 


T 


0.166 


0.17 


0.02 


0.15 


0.02 


0.15 


0.02 


0.22 


0.02 


0.22 


0.02 


0.2 


0.02 


A s 


0.9 


0.3 


0.04 


0.28 


0.04 


0.27 


0.04 


0.4 


0.04 


0.39 


0.04 


0.35 


0.04 



TABLE IV: Long term predictions assuming Planck for CMB, SNAP+NSNF for SNe and Planck/LSST for the CMB/LSS cross- 
correlation (X) for the "conservative" and "goal" cases. The SNe and X results are shown with priors from CMB temperature 
only (T) and CMB temperature and polarization data combined (TP). 







Model I, la errors 


Model II, la errors 






CMB 


SNe 


X(conserv.) 


X(g 


;oal) 


CMB 


SNe 


X(conserv.) 


X(g 


;oal) 


p 


fiducial 
v».ln.flfi 


T 


TP 


T 


TP 


T 


TP 


T 


TP 


T 


TP 


T 


TP 


T 


TP 


T 


TP 


Wo 


-0.98 


0.86 


0.46 


0.08 


0.08 


0.41 


0.33 


0.26 


0.24 


















Sw 


-0.29 


2.8 


1.5 


0.28 


0.27 


1.3 


1.0 


0.84 


0.76 


















(w) 


-0.9 


















0.14 


0.04 


0.03 


0.02 


0.04 


0.03 


0.03 


0.03 


A w 


-0.6 


















1.5 


0.39 


0.32 


0.24 


0.44 


0.3 


0.24 


0.19 


A z 


0.3 


















8.3 


4.5 


0.9 


0.79 


2.4 


2.2 


1.3 


1.2 


Q.M 


0.3 


0.01 


0.008 


0.007 


0.005 


0.01 


0.008 


0.01 


0.008 


0.01 


0.008 


0.009 


0.007 


0.01 


0.008 


0.01 


0.008 


h 


0.69 


0.01 


0.009 


0.006 


0.005 


0.01 


0.009 


0.01 


0.009 


0.01 


0.009 


0.01 


0.008 


0.01 


0.009 


0.01 


0.009 


10 3 u h 


24 


0.21 


0.16 


0.2 


0.16 


0.21 


0.16 


0.21 


0.16 


0.22 


0.16 


0.2 


0.16 


0.21 


0.16 


0.21 


0.16 


n s 


0.99 


0.005 0.004 


0.005 


0.004 


0.005 


0.004 


0.005 0.004 


0.005 


0.004 


0.005 


0.004 


0.005 


0.004 


0.005 


0.004 


T 


0.166 


0.04 


0.006 


0.03 


0.005 


0.03 


0.006 


0.03 


0.006 


0.05 


0.007 


0.03 


0.006 


0.03 


0.006 


0.03 


0.006 


As 


0.9 


0.07 


0.01 


0.06 


0.01 


0.06 


0.01 


0.06 


0.01 


0.1 


0.01 


0.06 


0.01 


0.06 


0.01 


0.06 


0.01 



in the absence of any priors can be found in Table HH 
Planck and SNAP are by far the most dominant infor- 
mation sources for these models, with cross-correlation 
adding a relatively weak contribution. 

For Model II we start by presenting the individual er- 
ror contours for the three experiments (Fig. I15f) . with 
the purpose of explicitely showing degeneracy directions. 
We will see that, once we consider its combination with 
Planck, the 

relative utility of the cross-correlation is significantly 
higher for Model II. The contours in Fig. El show that 
the cross-correlation can provide constraints on (w) and 
A w that are competitive with those from SNAP. The 
improvement on Planck is especially strong without the 
polarization data. The full list of parameter constraints 
is given in Table Hvl Our results confirm that SNe mea- 



surements are more suitable for constraining variations 
in w(z) at low redshifts. At higher redshifts, the cross- 
correlation becomes useful and can provide an indepen- 
dent constraint A„,. 

Figure El illustrates how a different choice of Model 
II fiducial parameters affects the constraining capabili- 
ties of the different probes. The alternative model has 
A w = —0.2 instead of —0.6 3 . In this case we note the 
SNAP+Planck constraint for A w is noticeably tighter 
than that from the Planck/LSST correlation. Also note 



3 The choice of the fiducial value of (w) was not modified as the 
transition time would be pushed to very high redshifts if it ap- 
proached -1. On the other hand, (w) smaller than —0.9 are 
strongly disfavored by current data. 
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FIG. 14: The expected la projected contours for the Model I 
parameters from Planck (thin blue), and Planck information 
combined with the Planck/LSST(goal) correlation (thick solid 
red) and SNAP supernovae (black dot). (Compare this to 
Fig, mi . 



FIG. 16: The expected la projected contours for the Model II 
parameters from Planck (thin blue), and Planck information 
combined with the Planck/LSST(goal) correlation (thick solid 
red) and SNAP supernovae (black dot). (Compare this to 

Fig.EU. 




FIG. 15: The expected la projected contours for the Model II 
parameters from the three probes: Planck (thin blue), 
Planck/LSST(goal) correlation (thick solid red) and SNAP 
supernovae (black dot). 



that this is mainly due to improvement in the sensitivity 
of the SNe, as the transition in this model happens at 
a lower redshfit. The size of the error on A w from the 
cross-correlation is only marginally increased compared 
to Fig. US] 



V. CONCLUSIONS 

We have studied the potential of the CMB/LSS corre- 
lation for constraining the evolution of the DE equation 
of state. We proposed a new parameterization, a "Kink" 
model which has the average equation of state (w) as an 
explicit parameter. A direct control on the value of (w) 
preserves the tight constraints on w obtained from the 
distance to last scattering measurements, while allowing 
w(z) to vary in other ways. 

Using Fisher analysis, we have made forecasts of the 
expected uncertainties in DE parameters extracted from 
the upcoming CMB, SNe and CMB/LSS correlation 
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data. We have considered both short- (WMAP, DES, 
SNLS) and long- (Planck, LSST, SNAP) term prospects. 
We find that in long term the cross-correlation will pro- 
vide competitive constraints on the evolution of w(z), 
in particular, on the total change in w, provided a suf- 
ficiently generic parameterization is used. We also note 
that the transition length parameter, A z , will not be well 
constrained even with the long-term combined data. 

Our results encourage further investigation of the 
CMB/LSS correlation as a probe of dark energy. The 
cross-correlation has a clear weakness - a large variance 
that hinders its utility as a precision cosmology probe. 
However, it also has certain attractive features, that in 
some cases can out weigh the lack of accuracy. These 
features are: 

• The cross-correlation probes the rate of the growth 
of structure as well as the growth. 

• It is practically independent of reionization details 
(r) and the details of the initial spectrum (n s ). 

• It probes the large scales (0.01 < k < 0.1/iMpc -1 ), 
where the evolution of structure is safely inside the 
linear regime and there is no need for higher order 
corrections. 

• On linear scales contributing to the ISW/LSS cor- 
relation the galaxy bias is expected to be indepen- 
dent of scale. 

• Cross-correlation is linear in bias, compared to the 
matter power spectrum that is quadratic. In fact, 
if one used a cross-correlation estimator that is 
normalized to the auto-correlations, that estima- 
tor would be independent of bias and the overall 
normalization of the primordial spectrum. 

• Cross-correlation can probe the evolution of dark 
energy at relatively high redshifts, inaccessible by 
luminosity distance measurements. 

In this work we have restricted our attention to 
Quintessence models for which the effective sound speed 
is of order of unity, so that the dark energy perturbations 
only affect the evolution of the largest structures. How- 
ever, if the dark energy clusters on smaller scales, it may 
leave a distinctive signature which can be detected with 



cross-correlation measurements. Although current data 
do not provide any constraints on the dark energy sound 
speed |59l l60j , the constraints from the next generation 
of galaxy surveys have been studied in |5fj| . 

We have only considered flat FRW models. Without 
the flatness assumption, CMB constraints on fijvf and 
h weaken because of a well-known degeneracy also 
relaxing constraints on w(z). We will explore constraints 
on DE using ISW without the flatness prior in a future 
work. 

The Model II type parameterization of w{z) consid- 
ered in the paper is not necessarily the most optimal for 
underscoring the power of cross-correlation constraints 
(although clearly better suited for this purposed than 
Model I). A more comprehensive study using a princi- 
pal component approach l63j | is a subject of ongoing 
research. 

Another potential use of the ISW cross-correlation is 
to measure the amplitude of the scalar perturbations and 
indirectly constrain the amplitude of the tensor modes 
when combined with the normalization of CMB spectra 
[64| . This is because the distribution of the large scale 
structures is not correlated with a stochastic background 
of gravitational waves. It has been also suggested j6^ 
that CMB/LSS studies may help to differentiate between 
Quintessence and models with modified gravity. 

For a fixed total number of galaxies the information 
extracted from the ISW studies depends on the choice 
of the redshift binning. The role of cross-correlation in 
constraining w(z) may be further increased by optimiz- 
ing the bin selection. Such an optimization is a topic of 
ongoing research. 
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